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Abstract 

The dependence of the dressed quark propagator on the quark chemical 
potential is investigated in various models based on the Dyson-Schwinger 
equations. We find that the critical chemical potential of the chiral symmetry 
restoration transitions is strongly dependent on the nature of the interaction 
kernel in the infrared region. 

The study of the thermodynamic phase structure of QCD has many implications for 
astrophysics and relativistic heavy ion collisions. An understanding of the behaviour of 
strongly interacting matter in regions of non-zero temperature and density will give insight 
into the nature of primordial baryogenesis and the dynamics of neutron stars. Heavy ion 
colliders such as RHIC and LHC will soon provide experimental data in thermodynamic 
regions where chiral symmetry is predicted to be restored. 

The structure of QCD at non-zero temperature and zero density has been explored ex- 
tensively using numerical simulations of lattice QCD. These results and other model based 
calculations suggest that there is a critical temperature, T crit ~ 170-190 MeV 0, above 
which normal nuclear matter undergoes a phase transition to a quark-gluon plasma in which 
chiral symmetry is restored. Unfortunately, current lattice simulation methods have signifi- 
cant difficulties when a non-zero quark or baryon density is introduced 0. For this reason, 
studies in models that implement the pertinent features of QCD are important. Here, we 
present results from such a study utilising the Dyson-Schwinger equation (DSE) framework 

At present we make no attempt to incorporate the colour superconducting ground state 
that has been shown to occur at asymptotic quark densities j| . The relevance of these results 
to QCD at densities of interest in nuclear physics (where perturbative results are suspect) is 
not well understood. These densities have only been investigated in simple approximations 
such as the Nambu-Jona-Lasinio model || and the Munczek-Nemirovsky (DSE) model || 
which may not be reliable representations of QCD in this regime. An extension of the current 
work may address these deficiencies and is under investigation. 

In principle, the DSE's determine all possible information about a quantum field theory 
by providing all of its Schwinger functions. Unfortunately, they form a countably infinite set 
of coupled integral equations with the equation for an n-point Schwinger function depending 
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on {n + 1) and higher point functions. Our approach is derived from these equations by 
specifying a scheme for truncating this infinite system through modelling the requisite, 
undetermined Schwinger functions. 

Previously, simple DSE models [|7||| have been used to investigate the restoration of 
chiral symmetry as the quark density (or chemical potential) increases. A phase transition 
has been found at a critical chemical potential, /i cr i t ~ 300-400 MeV. Here we examine the 
effect on this result of using different quark self-energy kernels. 

In Euclidean space (with metric g^ u = 5^, and Clifford algebra {7^,7^} = 25^, 
75 = —71727374), we separate the inverse quark propagator, S, at non-zero quark chem- 
ical potential, /i, into its spacelike and timelike vector parts, A and C, and scalar part, 
B: 

S-\p; = ij ■ pA(P; + ilAVpCfa + B(p- C) (1) 

where the quark momentum p = (p, uo p = p A + ifi). The propagator is renormalised at the 
(0(4) invariant) scale ( which is large enough that 

s-\p; 7° 2 = 5S" x (p; = *7 • p + m R (C) , (2) 

where Sq is the free propagator for a quark with running current mass mji((). In terms of 
the regularised self-energy E'(p, A), the inverse propagator can be expressed as 

S-^p; C) = Zfij ■ p + Z 2 (z 1a uj p + m bare (A)) + £'(p, A) , (3) 

where mbare(A) is the bare, regularisation dependent mass occuring in the Lagrangian, A a 
regularisation parameter (we use J* A = J A J dil 3 with a simple cutoff regulator A), and 

Z2 and Z 2 are respectively, the spacelike and timelike field renormalisations, dependent on 
both the regularisation and renormalisation scales. As we renormalise at \x = where 0(4) 
invariance is preserved, Z 2 = Z 2 . 

The DSE satisfied by the self-energy is 

E'(p, A) = g 2 [ A Df v {p - k)t a ^S(k)Tl(p, k) , (4) 

where D® b u (q) and T b v {p,k) are the full, nonperturbative gluon propagator and the quark 
gluon vertex respectively and {a,b} are colour indices with t a = 4r for the standard Gell- 
Mann SU(3) representation A a Q. 

As is the case for the propagator, the self-energy can be expressed in terms of three scalar 
functions Y/ B and E' c as 

S'(p; A) = Z7 ■ pS' A (p; A) + ij 4 0J p E' c (p; A) + S' B (p; A). (5) 

These functions satisfy the coupled DSEs, 

Z' x (p; A) = f A g"Df u (p-~k) itr [v xl ,t a S (k)T b u {p M , (6) 

where X = A, B, C; Va = ~{Z\V) • p/|p| 2 ), Vb = Z±, Vc = — (^1^74/^p) and the trace 
is over Dirac and colour space. We constrain the vertex renormalisation Z\ to be the same 



2 



as the field renormalisation, Z\ = Z2, which is consistent with the use of the bare vertex 
The renormalisation condition (^j implies that 



Z 2 « 2 ,A 2 ) 
m R (C) 



1-^(C,A)U 
^(C 2 ,A 2 )m bare (A) + S^(C,A)| / , 

Thus the regulator independent propagator functions are defined by 
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• 



A(p;C) = l + E^(p;A) 
B(p;Q=m R (0 + Z' B {p;A) 



(p;A)|g 



,2— £2 

b(p;A)|^ 



(8) 



C(p-Q = l + V c (p- A)-£' c (p;A)| 



=0 

p 2=C 2 



Equations ([I]) through (§) form a coupled set of nonlinear integral equations for the quark 
propagator, and once the self-energy kernel DY is specified we can proceed to solve them. 

Here, we consider various models for this kernel, and in particular, separately for D°^ v and 
. There are some constraints on the models that we can use to guide our choice; the main 
requirement is that they must be capable of producing dynamical chiral symmetry breaking 
(D%SB) and confinement || as this is what is found experimentally. The behaviour of these 
objects in the ultraviolet region is also dictated by the results of perturbation theory. 

The full Euclidean space gluon propagator can be written in Landau gauge as 



g 2 D; b M) 



g 2 S ab DM = 5 



QnQu\ Q{q 2 



Q 



Q 
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where we have assumed that the effect of quark chemical potential on the gluon propagator 
through quark loops is small in comparison to that on the quark propagator. To provide a 
model for it, we specify As the first example, we choose: 



G{q 2 ) 



(27r) 4 £5 4 (g), 



(10) 



which corresponds to the Munczek-Nemirovsky M (MN) model. This model has been studied 
in the context of finite temperature and density ||[10| and has the advantage that the integral 
equations of the DSE reduce to algebraic equations (albeit nonlinear). 
The remaining models are variations on the form: 



Q{q 2 



87r 4 A^ 4 (g) +4tt 2 (2- A) 



Gq 2 



LL- 



9 2 



47r 2 7 m 






g2(l)ln 


V Iy QCD J 


2" 

1 



:n) 



where the one loop anomalous dimension 7 m = 11N 12 2N , the QCD scale parameter is set 

c f 

at Aqcd = 0.275 GeV, and r = e 2 — 1. The last term in Eq. (|TTD implements the results of 
perturbative QCD to one loop. However, the quantities we study in this paper are deter- 
mined primarily by the low momentum structure of the kernel and are relatively insensitive 
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to this ultra-violet behaviour. The first two terms of the model provide some infrared en- 
hancement to the self-energy, leading to dynamical chiral symmetry breaking (D%SB). The 
exact form of the IR enhancement of the quark self-energy kernel is unknown and by 
varying A and to some extent u it is possible to explore various possibilities. For A = 2, 



the propagator QTTj) essentially reduces to that considered in Refs. [0,|T3j, however we have 
also included the UV tail which provides the correct leading-log asymptotic behaviour for 
the quark propagator. For A = 1 this is the propagator considered at zero density in Ref. 
|i~4| . Finally, with A = we get the propagator used in Ref. 

In order to motivate the different models used here for the vertex, we first briefly outline 
the method used to obtain the critical chemical potential. For a given self-energy kernel, the 
DSE's for the quark propagator are solved self-consistently to find solutions representing the 
Nambu-Goldstone phase (characterised by D%SB) and a Wigner-Weyl phase (corresponding 
to the quark- gluon plasma where chiral symmetry is restored). The energetic stability of 
these two solutions is then compared using an effective action for the composite operators 
S and D. 

One commonly used truncation of the DSE (|6|), called the rainbow approximation, in- 
volves replacing the full vertex in Eq. (|^) by the bare vertex, 

K(KP)=ta7u- (12) 
In this case, the correct action to use is the Cornwall- Jackiw-Tomboulis (CJT) effective 



action [16[] which is given by 

AcMS] = TrLn [S^S] + ±Tr [ES] . (13) 

We wish to explore the effects of moving beyond the rainbow approximation (ideally to 
a vertex that respects the various symmetries of the full vertex, see Refs. @,0]), but must 
first address the question of which effective action should be used to measure the energy of 
the different solutions. For a more complicated vertex, functionally dependent on the quark 
propagator, the CJT action is not the correct object to use [TB|. However the correction 
required is so far only calculable in a limited number of casesQ. This restricts the types of 
models that we can use for the vertex. 

Consequently, we use the next level of truncation of the DSE (|6]); a one loop vertex 
involving the non-perturbative quark and gluon propagators. This vertex was introduced 
in Ref. |B| and used further in Refs. P,|2"0[|. In combination with the appropriate Bethe- 



Salpeter kernel, it preserves the axial-vector Ward identity and is systematically improvable. 
The form for this vertex is 



1 Indeed, recent lattice calculations [|ll]] and DSE studies [12] suggest the enhancement may in 
fact not occur in the gluon propagator but through the vertex. In our approach, it is full kernel 
DT that is modelled and the particular separation into gluon propagator and vertex function is 
arbitrary. 

2 This is essentially the same problem that frustrates a consistent truncation of the DSE and 
corresponding BSE equations 
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The correct action to use in combination with this vertex is, 

A 1L [S,D]=A c MS}+A r [S,D], 



where [18], 



Ar[3, D\ = -kj d A pd 4 qd A rD pa (p - q)D^{r - g)Tr [^S^Srf S rl ° S, 



p-q+ri 



(14) 



(15) 



(16) 



Having described the basis of the calculations, we now turn to the most important results. 
For the MN model with a bare vertex the analytic solutions are known. In the chiral limit, 
the two relevant solutions are: 



Si : 



B(p) = 



A(p) 




S 2 : B(p) = 2VG=P Z A(p) = 2. 



(17) 



In order to define the symmetric and D^SB solutions for the MN model at non-zero quark 
chemical potential, it is necessary to provide a prescription for the transition from solution 
S 2 to solution Si as \p\ — > oo. To generalise the \i = case (where the reality of B determines 
the transition point), we set <Se>xSB = S 2 for 9ftQo 2 ) < G + /i 2 and Sd^sb = <Si otherwise. We 
stress that while this prescription is a reasonable generalisation of the model to finite fi, it 
is essentially arbitrary and different results are obtained with other prescriptions fll8"| . 

For this model, the critical chemical potential (/i cr it such that Acjt = ^■cjt[5'dxSb] — 
-4cjT[>S'symm.] = 0) is displayed as a function of the coupling strength Q in Fig. 1. As Q is 
increased, there is a corresponding monotonic increase in /i crit which can be fit exceedingly 
well with the simple form, /i cr it({?) = 0.335v^- The "physical" value of the model parameter 
Q can be set by requiring that it yield the correct zero density pion decay constant, /„-. For 
this coupling, the critical chemical potential is /i cr i t = 340 MeV. 

For the more phenomenologically appropriate models using the Maris-Roberts type ef- 
fective quark- ant iquark interactions, Eq. (|TTD , the integral equations @ are solved using an 
iterative procedure for ( = 19 GeV and A ~ 10 3 GeV. The quark condensate is extracted 
from the asymptotic behaviour of the quark mass function M(p 2 ) = B(p 2 ; ()/A(p 2 ; (), 



M(p 



2\ P 



I In 



A 2 

QCD 



1-7" 



(18) 



which is accurate up to 0( 



A 2 

^QCD ' 

c 2 



The zero density pion mass and decay constant are cal- 



culated using analytic approximations to the solutions of the corresponding Bethe-Salpeter 
equation |Tj|. These values are used to fix the bare quark mass and the parameters Q and 
lj in the gluon propagator. 

Calculations of the critical chemical potential in these models give qualitatively similar 
results to those of the MN model. Data from the A = 1 and A = 2 cases is shown for 
a range of Q in the upper panel of Figure 1. It is interesting to explore the parameter 
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dependence of these models. Figures 2 and 3 present the normalised CJT action differences 
— A cjt [/^]/Acjt[0] as a function of \i for variation of the three parameters Q, A and uj. 
The results for different parameter choices are normalised to —1 at zero chemical potential 
as |Acjt[0]| increases rapidly with Q. For A = 2, Fig. 2 shows that as the total infrared 
strength Q is increased yn C rit increases in a manner similar to that of the MN model. Fig. 
3 shows the normalised action differences for variation with A (left panel, Q = 0.26 GeV 2 , 
uj = 0.3 GeV) and uj (right panel, Q = 0.26 GeV 2 , A = 1). As A decreases from 2 
(where the entire infrared strength is in the delta function at zero momentum) to (where 
the effective quark- ant iquark interaction vanishes at zero momentum), the critical chemical 
potential decreases for a given coupling strength and fixed value of uj. Similarly, as the 
infrared strength of the propagator is broadened by increasing uj for a given Q and A, /i crit 
again decreases. Table 1 summarises the parameters used and the corresponding values of 
the critical chemical potential for the various models considered here. 

The effect on fi crit of using the one loop vertex is twofold. First, the solutions themselves 
are modified, resulting in a change to both the physical values for the model parameters 
and the action differences. Using the corrected action, Ail, and the MN gluon propagator, 
the critical chemical potential is f^ t (G) = 0.279\/^. Secondly, the correction to the action 
produces an additional shift. We illustrate this for the MN model though similar, (but 



less clean) conclusions can be drawn from the Maris- Roberts type models [18|. Figure 1 
(dashed line) shows the curve of critical chemical potential for the one loop vertex MN 
model evaluated with the full action, Ail- Also illustrated in the bottom panel of this figure 
is the effect of the correction term, Ar, on the critical chemical potential. Data are shown 
for the difference between the critical chemical potentials obtained from the one loop vertex 
solutions using the CJT and full, one loop actions; specifically for 5/x crit = [/4rit T — A*crit] /Vcrit- 
It can be seen that the effect of the correction term is insignificant (|<5/i cr it| < 2% for the 
range of values for Q considered here) in comparison to the modification of /j, crit because of 
the use of one loop solutions. 

In summary, the studies performed here make it clear that the critical chemical potential 
is an "observable" that is sensitive to the non-perturbative structure of the quark self energy 
kernel. For reasonable models with parameters fitted to the quark condensate, pion mass and 
decay constant, we find that /Vit lies between 300 and 650 MeV. Our studies utilising a one 
loop vertex support the possibility that the CJT action, although not strictly appropriate, 
provides a reliable extraction of // cr it for vertices other than the bare vertex. Whether this 
remains true for (physically) more acceptable vertices such as those of Ball and Chiu or 
Curtis and Pennington |L7| remains an open question. 
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TABLES 





G \C1pV 2 ] 


UJ 1 UC V J 


m \A/Tt>V\ 
I itq yivi c v J 


\ \1Q)) [o-ev j 




/ / c-Tj- 1 IVlC V J 


Merit L lvi c V J 


MN(bare) 


1.0 




17 


0.216 


93 


138 


340 


MN(bare) 


1.20 




17 


0.236 


102 


145 


370 


MN(loop) 


0.93 




15 


0.211 


93 


137 


270 


MN(loop) 


1.15 




15 


0.236 


105 


144 


300 


MR(A = 2) 


0.82 




9.8 


0.236 


90 


138 


630 


MR (A = 1) 


0.61 


0.3 


6.6 


0.236 


77 


138 


610 


MR (A = 0) 


0.45 


0.3 


5.7 


0.236 


75 


139 


530 



TABLE I. Parameters and observables are shown for a variety of models studied. The quark 
masses, m q , listed for the Munczek-Nemirovsky (MN) models are the bare masses as this model 
is not renormalised. For the Maris-Roberts (MR) type models the value listed is the renormalised 
mass at 19 GeV 2 . The different cases for each particular model are when either the condensate or 
pion decay constant is fitted. 
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FIGURES 
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FIG. 1. In the upper panel, the calculated critical chemical potential is shown for the various 
self-energy kernels as a function of the coupling strength, Q (G in axes labels). The cases are: (open 
squares) Munczek-Nemirovsky gluon propagator and bare vertex (CJT action); (open triangles) 
MN gluon and one- loop vertex (calculated using one- loop action); (open circles) Maris- Roberts 
(MR) type propagator, bare vertex with A = 1 and co = 0.3 GeV; (open stars) MR propagator, 
bare vertex with A = 2. Simple fits of the form // C rit(£7) = a are shown for each data set. In the 
MN cases, the best fits occur with (3 = 1/2 and with a = 0.335(0.279) for the bare (loop) vertex. 
With the MR propagator the best fits are: (A = 2 case) a = 0.696, (3 = 0.516 and (A = 1 case) 
a = 0.849, (3 = 0.656. In the lower panel, the effect on the critical chemical potential of using 
the corrected action for the one loop MN model solutions is shown. This effect is small and the 
magnitude of Sfi cr it = [^criF ~ A*crit]/ ^crit i s ^ ess than 2 % over the range of Q studied. 



9 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 

yu [GeV] 

FIG. 2. Variation of the normalised CJT action difference, — Acjt(a*)/Acjt(0), with chemical 
potential is shown for a range of Q (labelled in units of GeV 2 ). The critical chemical potential for 
a given Q is reached at the zero intercept of that curve. 




0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 

fx [GeV] \jl [GeV] 

FIG. 3. Variation of the normalised CJT action difference, — Acjt(m)/Acjt(0), with chemical 

potential is shown for a range of A (left panel, with fixed Q = 0.26 GeV 2 and u = 0.3 GeV) and 

oo (right panel, with A = 1 and fixed Q = 0.26 GeV 2 ). 
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